Inclusive Hadron Production in p+p Collisions at 

STAR 



Frank Simon (for the STAR collaboration) 

Massachusetts Institute of Technology 

Abstract. The STAR experiment at RHIC has measured a variety of inclusive hadron cross sections 
in p + p collisions at y/s = 200 GeV. Measurements of the differential cross section for inclusive 
charged pion production at mid rapidity and for inclusive neutral pion production at forward rapidity 
(3.0 < tj < 4.2) as well as the first preliminary result from STAR for the differential cross section for 
inclusive neutral pion production near mid rapidity are presented. These cross sections are compared 
to next-to-leading order perturbative QCD calculations and can provide constraints on the pion 
fragmentation functions. Good agreement between data and pQCD has been found for all three 
cross sections. 



INTRODUCTION 

One goal of the polarized p+p program at the Relativistic Heavy Ion Collider (RHIC) is 
the study of the gluon polarization AG in the proton via spin asymmetry measurements 
in a variety of processes [1]. Such studies typically involve fits to polarized data using 
perturbative QCD (pQCD) evaluation for the various processes. In order to validate the 
application of pQCD a good understanding of the (unpolarized) cross sections of these 
processes is crucial. In addition the study of inclusive cross sections of identified parti- 
cles provides important information on the general applicability of pQCD calculations 
in the RHIC energy regime as well as on fragmentation functions. 

The STAR experiment at the RHIC §] is well suited for inclusive hadron measure- 
ments over a wide range in rapidity due to its central large acceptance time projection 
chamber (TPC) and electromagnetic calorimeters, a time of flight (TOF) detection sys- 
tem for particle identification and the forward neutral pion detectors (FPD). 



EXPERIMENT AND DATA ANALYSIS 

The forward neutral pion cross section is measured with the FPDs These detec- 

tors are high resolution electromagnetic calorimeters consisting of lead glass crystals, 
covering the pseudorapidity range ~ 3.0 < \ f] \ <~ 4.2 on both sides of the interaction 
region. Events were triggered by a minimum energy requirement in the FPD in addition 
to the p + p minimum bias trigger based on a coincidence of the beam beam counters 
(BBC). Neutral pions are reconstructed by reconstructing the invariant mass from their 
two decay photons. The mass peak position is used to calibrate the FPD energy scale to 
a precision of ~1%. The data presented here were taken in the 2002 and the 2003 RHIC 




FIGURE 1. Inclusive hadron cross sections in p+p collisions at ^/s = 200 GeV. 

Left: Mid rapidity n° cross section: a) compared to NLO pQCD calculations with two sets of frag- 
mentation functions; b) relative statistical error (bars) and point-to-point systematics (band); c) relative 
difference to NLO pQCD calculations. For the KKP fragmentation function, the scale uncertainty is 
indicated by the dotted (jJ. = 2 pj) and dashed-dotted (fi = Pt/2) lines. 

Right: d) Forward rapidity n° cross section for three rapidity intervals [4], compared to NLO pQCD 
calculations, e) Mid rapidity charged pion cross section [6], compared to NLO pQCD with two sets of 
fragmentation functions, f) Ratio with theory for the KKP fragmentation function. 



p + p run, with an integrated luminosity of 350 nb . The cross section is determined in 
three ranges with mean pseudorapidities (rj) = 3.3, 3.8 and 4.0. 

Charged pions are identified using the TPC and the TOF system. The TPC tracks 
charged particles in the pseudorapidity range — 1 < 77 < 1 . For the 2003 RHIC run a first 
prototype of the TOF was installed and is used to provide pion-kaon separation below 
p t = 2.5 GeV [5]. Above that threshold, the ionization energy loss dE/dx in the TPC is 
used to identify pions. Due to the relativistic rise of the energy loss, the pion dE/dx is 
~ 15% higher than that of kaons and can be used to separate the two particle species. A 
detailed description of the charged pion analysis is given in [Jg]. 

Neutral pions near mid-rapidity are detected with the barrel electromagnetic calorime- 
ter BEMC [7]. For the 2005 run period, half of the BEMC was fully installed and com- 
missioned, giving a coverage of < r\ < 1 for all azimuthal angles <p. The BEMC is a 
lead scintillator sampling calorimeter with a granularity of 0.05 x 0.05 in Arj x A(p. For 
the reconstruction of high p t neutral pions the identification of the two decay photons 
with a small opening angle is crucial. This is achieved using the shower maximum de- 
tector SMD, a wire proportional counter with cathode strip readout at ~ 5Xq depth in 
the calorimeter modules which provides a segmentation of 0.007 x 0.007 in Ari x A<p. 



Neutral pions are identified from the 77 invariant mass spectrum. The mass and p t are 
determined from the opening angle between the two photons measured with the SMD 
and from the neutral energy deposited in the calorimeter. A charged particle veto re- 
jects calorimeter hits that have contributions from charged particles. Dedicated online 
triggers, selecting events with high energy deposit in a BEMC tower in addition to the 
minimum bias trigger condition, were used to select events containing high p t neutral 
pions. These high tower (HT) triggers used thresholds of ~2.8 GeV and ~3.6 GeV to 
allow prescaling of the abundant low-pT signals and thus improve the detector live-time. 
The absolute energy scale of the BEMC is calibrated in situ using high energy electrons 
whose momentum is determined by the TPC, while a relative tower by tower calibration 
is performed with more abundant minimum ionizing particles. For the preliminary n° 
cross section presented here, a subset of the available 2005 data with an integrated lu- 
minosity of ~ 44 /lb ' _1 for minimum bias and ~ 0.4pb~ ! for HT triggers was analyzed. 
The cross section is given by 
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where Ap t and Ay are the bin widths in p t and rapidity, Jzf is the sampled luminosity, 
Njc is the number of reconstructed n° in the bin and c is an overall correction factor 
that contains the corrections for acceptance, reconstruction and trigger efficiency thus 
giving the true number of pions from the number of reconstructed pions in a given bin. 
These corrections are determined from a full GEANT simulation of the detector, using 
p + p events generated with the PYTHIA event generator. Good agreement between the 
observed and the simulated 77 spectra has been found. 



RESULTS AND DISCUSSION 

Figure[T]shows the inclusive cross sections for neutral pions at mid- and forward rapidity 
and for charged pions at mid-rapidity. All cross sections are compared to NLO pQCD 
calculations using the Kretzer H] and the KKP ifioll fragmentation function sets. 
The pQCD calculations are performed with a factorization scale \x = p t , unless stated 
otherwise. 

Figure Qt) shows the preliminary cross section for inclusive 7T° production in the 
rapidity interval 0.1 < y < 0.9 together with the NLO pQCD calculations. Figure [TJ)) 
shows the relative statistical errors of the data points as error bars and the preliminary 
point-to-point systematic errors as the shaded band. The point-to-point systematics are 
dominated by errors on the yield extraction from the 77 invariant mass spectra. Figure 
Ut) shows the relative difference of the measured cross section to the NLO pQCD cal- 
culations for both the Kretzer and the KKP fragmentation functions. The shaded band 
shows the preliminary total systematic error, which is dominated by the energy scale 
uncertainty of the calorimeter calibration, estimated to be ~5% for the present data. The 
NLO calculations using the KKP fragmentation functions show excellent agreement 
with the data to p t below 3 GeV/c. The Kretzer fragmentation functions consistently 
underpredict the data. To indicate the scale uncertainty of the pQCD calculations, calcu- 



lations using the KKP fragmentation functions with scales of ji = 2p t and ji = p t /2 are 
shown by the dotted and the dashed-dotted lines, respectively. 

FigureUJi) shows the inclusive neutral pion cross section at forward rapidities in three 
pseudo-rapidity bins with (tj) = 3.3, (17) = 3.8 and (77) = 4.0. The data points are again 
compared to NLO pQCD calculations using the Kretzer and the KKP fragmentation 
functions. At low rapidities, the KKP fragmentation functions show much better agree- 
ment with the data than the Kretzer set. At high rapidities and lower 7T° energy the data 
tends towards the Kretzer fragmentation functions. One consideration here is that the 
transverse momentum in the high rapidity bins is quite low and where pQCD at NLO 
might not be a good expansion. For r\ = 4.0 a K° energy of 40 GeV corresponds to p t ~ 
1.5 GeV/c, while at 77 = 3.3 an energy of 40 GeV corresponds to p t ~ 3.0 GeV/c. 

Figure [1^) shows the inclusive charged pion cross section at mid rapidity, compared 
to NLO pQCD using the two different sets of fragmentation functions discussed here. 
Again, the KKP set shows much better agreement with the data than the Kretzer frag- 
mentation functions. Figured) shows the ratio of the data to the NLO calculations using 
the KKP fragmentation functions. The excellent agreement to p t as low as 2 GeV/c is 
apparent. At very low p t , the data tends to lie below the predictions of the KKP frag- 
mentation and more towards the Kretzer set, as also observed for forward neutral pions. 
However, in the low p t regime the scale uncertainties of the pQCD calculations are very 
significant, as apparent from figureUt), demonstrating the limited applicability of pQCD 
calculations in this regime of low transverse momentum. 

SUMMARY AND OUTLOOK 

The presented inclusive hadron cross sections in p + p collisions at ^Js = 200 GeV 
show very good agreement with NLO pQCD calculations down to transverse momenta 
well below 3 GeV/c. This demonstrates the applicability of these calculations for the 
interpretation of measured spin asymmetries to extract the gluon polarization in the 
proton. All three cross section measurements presented here consistently favor the 
KKP fragmentation function set over the Kretzer set for p t > 2 GeV/c. However, scale 
uncertainties of the calculations are of comparable size to the difference between the 
fragmentation function sets. Continuing analysis of the 2005 p + p STAR data-set will 
lead to improved statistics and reduced systematic uncertainties in the inclusive n° cross 
section at mid-rapidity. 
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